DFT+U study of magnetic order in doped La2Cu04 crystals 
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This article presents the results of several magnetic phases of doped La2-xSrxCu04 using density- 
functional theory with an added Hubbard term (DFT+U). Doping factors from x = to 0.25 were 
examined. We found that a bond centered stripe is the magnetic ground state for x = | and x = |. 
No stable stripe order was found for x = |. Analysis of the electron density revealed that apical 
oxygen atoms, those located above and below the copper atoms in the Cu02 planes, hold a non 
negligible part of the holes at large doping and present a small spin polarization. Finally, the charge 
reorganization caused by the magnetic stripe modulation was studied for bond centered and atom 
centered stripes. 



I. INTRODUCTION 

Over the last few years, studies on high-Tc cuprates 
have reported that the superconducting state competes 
with magnetic stripe phases For low doping, the su- 
perexchange^ is strong enough to favor an antiferromag- 
netic order. With increased doping, this interaction is 
weakened and a superconducting phase can arise. In the 
La2-a:Sra^Cu04 system, a sharp decrease in the super- 
conducting temperature is observed for a doping factor 
oi X = \. At that fraction of strontium substitution, 
static stripe ordered phases have been detected)^ Mag- 
netic phases were also identified for other dopings, up to 
0.18.-^ 




FIG. 1. Atomic structure of La2Cu04. Lanthanum atoms are 
represented by blue spheres, copper atoms by gold spheres 
and oxygen atoms by red spheres. Apical oxygens are located 
above and below the copper. In plane oxygens d and O2 are 
in X and y directions in the Cu02 plane. 

The atomic structure of the La2Cu04 is depicted in 
Fig. [1] The oxygen atoms in the middle of the octahe- 



dra and the copper atoms form the copper-oxide plans 
where stripe ordered phases have been reported. In or- 
der to understand the phenomena that yields spin stripes 
in La2-a?Sra;Cu04, two commensurate magnetic arrange- 
ments were proposed in previous works. ^ ^ They are de- 
picted in Fig. |2] and [3) Their types of magnetic structure 
break the translation symmetry of the unit cell. The 
main difference between these stripes is the location of 
holes in the structure. In the first one, the holes are 
located on copper atoms and it is generally referred as 
an atom-centered stripe (AC). The other configuration is 
known as a bond-centered stripe (BC) because holes are 
localized on oxygen atoms between copper atoms in the 
Cu02 planes. 

It is relatively difficult experimentally to measure the 
spin configuration in the stripe configurations. To help 
characterize these phases, we conducted an ah initio 
study on different magnetic structures to determine their 
relative stability. To this end, we used an extension of 
density- functional theory that include a Hubbard term 
within the hamiltoniaiP^ which we found is necessary 
to stabilize the stripe configurations. We also looked at 
these magnetic phases for different doping factors. Pre- 
cisely, we report the study of the stability of magnetic 
orders for doping ranging from x = to 0.25. Bond- 
centered and atom-centered stripes were considered as 
well as the antiferromagnetic state. The next section de- 
scribes the method used and it is followed by a section 
where we analyze the stability of the different magnetic 
stripes, the amplitude of the magnetic moments on cop- 
per atoms and the atomic charge density for each stripe. 
Finally, holes distribution is calculated in different stripes 
and compared to the antiferromagnetic phase with the 
same doping. It is found that for x=0.125 and 0.25, the 
ground state is the BC stripe and for large doping val- 
ues, a non negligible part of holes are localized on apical 
oxygen atoms. 



II. METHOD 

This study has been conducted within the Kohn-Sham 
approaclF^ of the density-functional iheoiy^^ (DFT) us- 
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FIG. 2. Representation of the bond-centered stripe. The oxygen atoms are represented by circles and the size of each circle 
is proportional to the hole content on this oxygen atom. For every copper atom, an arrow indicates the orientation of the 
magnetization of its 3(1^2 _y2 orbital. On top of the picture are presented respectively the copper magnetizations and the charge 
on Oi and O2. The values are reported with respect to their position in the stripe. The values were obtained using t/= 4 eV 
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FIG. 3. Representation of the atom-centered stripe. The oxygen atoms are represented by circles and the size of each circle 
is proportional to the hole content on this oxygen atom. For every copper atom, an arrow indicates the orientation of the 
magnetization of its 3d^2_y2 orbital. On top of the picture are presented respectively the copper magnetizations and the charge 
on Oi and O2. The values are reported with respect to their position in the stripe. The values were obtained using U= 4 eV 



ing the Projected Augmented Waves (PAWj^ formal- 
ism as implemented in the ABINIT code.'^^'^ For the 
exchange-correlation functional, we used the general gra- 
dient approximation of Perdew-Burkes-Ernzerhof.^ To 
account for the correlation energy of the strongly local- 
ized copper 3d orbitals, we employed a Hubbard term 
within the functional (DFT+U).^^ The Hubbard param- 
eters U and J within this method need to be assigned.'^ 
We have considered two values for the U term, 4 and 
8 eV, with a J term of respectively 0.4 and 0.8 eV. The 
value of 8 eV for U was chosen because it was determined 
using a density-functional approach^^ in a previous study 
of this system^ using a local density functional for the 
exchange-correlation energy. Since the present study 
uses a general gradient approximation for the exchange- 
correlation functional which is thought to account better 
for the correlation energy of localized orbitals, we also 
used the lesser value of 4 eV for the U term. A plane- 
wave basis set with a maximal kinetic energy of 35 Ha 
was employed. The Monkhorst-Pacli^ k-point grids for 
the 8 units supercell, 6 units supercell and 4 units su- 
percell were respectively 4x4x2, 6x6x4 and 8x8x6. Full 
relaxation of the structures was done in order to obtain 



their lowest energy configuration. 
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FIG. 4. DFT+U band structure of the EC stripe with 
U—4: eV. The band crossing the Fermi level corresponds to 
the direction perpendicular to the stripe. The other path is 
oriented along the stripe. 

The doping of the antiferromagnet La2Cu04 was done 
by subtracting a number of electrons from the system 
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and adding a neutralizing background. The atomic struc- 
tures studied were based on the tetragonal cell with space 
group I4/mmm. This structure gave good results for 
doped La2Cu04 crystals in a previous work.^^ It was 
also observed that for static stripe ordered phases with 
neodymium added to pin down the ^ niag netic order, the 
crystal shows I4/mmm symmetries'^ which validate 
the use of this atomic structure for the study of stripes in 
La2Cu04. Recently, it was observed that stripes can oc- 
cur in crystals doped only with barium.^ This indicates 
that magnetic rare-earth atoms are not required for the 
appearance of stripes in La2Cu04 crystals further vali- 
dating our doping procedure. 
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TABLE I. Energy differences per copper atom (in eV) between 
the magnetic phases considered and the normal metallic phase 
for three doping values. The positive values mean that the 
magnetic phases are more stable than the normal phase. Re- 
sults are given for two U values. The antiferromagnetic phase 
is named AF. Note that for x = |, only one type of stripe 
was calculated. See text for details. 



8 eV does not change the metallic property of the crys- 
tal but it reduces the bandwidth of the highest occupied 
band from 0.3 to 0.2 eV. It also increases the band gap 
between the half filled band and the first completely un- 
occupied band to 1.9 eV. This is similar to the non-doped 
antiferromagnetic insulator phase where an increase of U 
pushes the unoccupied band to higher energy. Further- 
more, we note that the atom-centered stripe is the second 
lowest energy state. Only a small energy difference sep- 
arates the BC stripe from the AC stripe for both values 
of U used. We also see that the antiferromagnetic state 
is barely stable as compared to the normal state for a 
U value of 4 eV but that the stability of this state is 
greatly increased for a /7 = 8 eV and it is comparable 
to the AC stripe total energy for this value of U . The 
reason is that an increase of U favors the magnetization 
of the electrons in the 3dj;2_y2 orbital. Finally, we notice 
that the difference in energy between of the BC stripe 
and the AC stripe increased when the U term changed 
from 4 to 8 eV. 



III. RESULTS AND DISCUSSION 

First, we will present the study of the stability of mag- 
netic phases and the copper magnetization. It will be 
followed by the charge analysis of the stripe phases. Fi- 
nally, the hole distribution is computed for the different 
stripes considered. 



A. Total energy stability analysis 

For X = 0.125, we constructed a 8 units supercell 
in which we removed one electron. We first calculated 
the normal state given by the DFT+U approach with- 
out spin polarization. This state is clearly metallic due 
to the odd number of electrons per primitive cell. We 
then computed the three magnetic states, the AC and 
BC stripes and the antiferromagnetic phase, and we cal- 
culated the total energy differences between these states 
and the normal metallic state. These results are reported 
in Table [Tj We find that the BC stripe has the lowest 
energy for both values of U used. This indicates that 
the bond-centered stripe should be the ground state of 
Lai.875Sro.i25Cu04 which is in agreement with previous 
results reported in the literature.^ The band structure of 
this stripe for U = 4 eV, shown in Fig. [4j shows a metal- 
lic behavior. The band crossing the Fermi level is in the 
direction perpendicular to the stripe. Increasing U to 
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FIG. 5. Representation of the two magnetic stripe orders. 
The upper panel represents the bond-centered stripe (BC) 
and the lower panel the atom-centered stripe (AC). The cir- 
cles represent the oxygen atoms. Arrows, located on the cop- 
per sites, indicate the orientation of the magnetization of the 



3d^ 



orbitals. 



Using the same approach, we studied other magnetic 
configurations with different doping levels. First, we con- 
structed a supercell containing 4 primitive cells. Remov- 
ing an electron from the total density, we obtained a dop- 
ing fraction per formula unit of |. We then computed 
three commensurate magnetic structures, one being the 
antiferromagnetic state. The other magnetic structures 
are shown in Fig. |5j These magnetic orders are based 
on the 8 units supercell spin configurations. They repre- 
sent possible states for Lai.75Sro.25Cu04. The difference 
in total energies for these different structures are also 
reported in Table [l| Again, the BC arrangement is the 
lowest energy state. Nonetheless, the energy differences 
are smaller compared to the case of x = h. This is due to 
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hole formation in oxygen and copper orbitals that Umits 
the effective spin coupling on copper atoms. Neverthe- 
less the DFT+U method predicts that a magnetic order, 
similar to the BC structure, can be found in La2Cu04 for 
a doping of x = 0.25. Interestingly, the AC stripe is not 
favored compared to the antiferromagnetic state for this 
doping and U =8 eV. For the smaller value of [/, the AC 
phase is more stable than the antiferromagnetic order. 
This is explained by the important increase in stabil- 
ity of the antiferromagnetic phase for U =8 eV compared 
to t/ = 4 eV. For this doping factor, the superconduct- 
ing temperature is almost zero. Consequently, a stripe 
phase would not generate an important drop in super- 
conducting critical temper at ure(Tc). Besides, the results 
obtained by the DFT+U show that for a = 4 eV in 
copper orbitals, the disparity between the normal metal- 
lic and the antiferromagnetic states is 2 meV per copper 
atom. Because no static antiferromagnetic state is ob- 
served for X = 0.25, the small difference in energy ob- 
tained with this value of U seems to describe correctly 
this system. The BC stripe being an order of magnitude 
more stable than the antiferromagnetic state, it means 
that BC stripe could be detected in Lai.75Sro.25Cu04. It 
is possible though that this state is fluctuating, like in the 
case of stripe order in 0.125 doped La2Cu04. The results 
with U =8 eV show that for a greater Coulomb repulsion, 
the stripe ordered phase is even more favored. However, 
for this value of doping, it is reasonable to think that 
the effective U term would be even smaller than 4 eV 
which would explain why the metallic state is the most 
energetically favorable state at that doping. 

In order to investigate the stability of magnetic order 
near the optimal doping {x ~ 0.19), we studied a 6 units 
long supercell with a doping factor of ^. This time, mag- 
netic orders based on the stripes for x = 0.125 were not 
possible without breaking the commensurability of the 
magnetic unit cell. As a result, we considered a magnetic 
stripe consisting of three up spins followed by three down 
spins. We found that this phase was not stable com- 
pared to the antiferromagnetic state for the two values 
of U. The AC and BC stripes being incompatible with 
the commensurability of the supercell, DFT+U predicts 
that no commensurate stripe order should be detected 
for X = ^. In this case, the energy difference between the 
antiferromagnetic state and the normal metallic state is a 
little higher than for the case X — 0.125. The difference 
in doping value being close to 0.125, the antiferromag- 
netic phase of X = I should be similar to 0.125 phase. 
This will be shown in the next section. 



B. Magnetic moments 

For each doping value, the main difference between 
the stripe configurations is the region where holes are lo- 
cated. Depending on where the holes are based in the 
Cu02 plane, it modifies the interactions between copper 
and oxygen atoms. As mentioned before, it affects the su- 



perexchange interaction resulting in a modulation of the 
magnetic moments throughout the supercell. To analyze 
the magnitude of the moments, we can split the stripes 
in antiferromagnetic and ferromagnetic parts. The BC 
stripe can be described by alternating antiferromagnetic 
and ferromagnetic regions. On the other hand, the AC 
stripe contains two antiferromagnetic domains separated 
by non-polarized walls. We found that, for x = 0.125, see 
Fig. [2] and [3) the antiferromagnetic parts for the BC as 
well as for the AC stripes have magnetic moments very 
close to the antiferromagnetic phase with the same dop- 
ing factor. In fact, the copper magnetization in the anti- 
ferromagnetic phase are 0.45 and 0.63 /i^ for = 4 eV 
and 8 eV respectively. At zero doping, the experimental 
value is 0.6 jub which is well reproduced with a /7 of 8 eV, 
however, for this value of the Coulomb repulsion, the an- 
tiferromagnetic doping region is too large. The variation 
in magnetic moments for the BC stripe is 0.12 /hb for 
/7 = 4 eV (0.05 /i^ for /7 = 8 eV). Nonetheless, each 
copper atom has a non-zero magnetization. For the AC 
stripe, this variation is greater because of the two non- 
polarized atoms. But the difference of absolute values 
for magnetic moments between polarized atoms in the 
AC stripe is very low. Finally, we see that the DFT+U 
approach can give rise to different modulations of spins in 
La2Cu04 doped crystal. For x = 0.25, the antiferromag- 
netic phase has a magnitude of spin polarization on the 
copper atom of 0.31 fiB for a Coulomb repulsion of 4 eV. 
This value is lower than those for x = 0.125. For ^ = |, 
we found a magnetization on the copper atom of 0.47 jub 
which is as high as for x = 0.125. This is consistent with 
the high stability of the antiferromagnetic phase of the 6 
units supercell. 



C. Charge analysis and hole distribution 

The spin configurations in the supercells can be un- 
derstood by the hole distribution, i.e. the localiza- 
tion of holes in the Cu02 planes. In order to analyze 
this distribution through the crystal, we calculated the 
Bader charges^^ in these systems with a of 4 eV. This 
method splits up the electronic density associated to a 
specific atom. Bader volumes, containing atomic do- 
mains, are defined by zero flux surfaces of the electronic 
density. First, we computed the antiferromagnetic phases 
for < X < 0.25. As mentioned before, with a value 
of 17=4 eV, the DFT+U can find an antiferromagnetic 
phase even for high doping {x = 0.25). For lower /7, 
the predicted antiferromagnetic doping region can be re- 
duced, however, the resulting moments for the non-doped 
system will be much lower than the experimentally ob- 
served value. It seems that within the DFT+U approach, 
the U value should be lowered as doping increases to 
account for the experimental phases observed. In the 
present study, we considered a U term of 4 eV for all 
doping, which seems a reasonable value. To obtain the 
location of holes in the crystal, we take the difference be- 
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tween the non-doped system and the doped one. Fig. |6] 
shows the increase of hole content on copper, apical oxy- 
gen, in-plane oxygen and lanthanum atoms with respect 
to the doping factor x. When adding up the hole fraction 
of each atom in the primitive cell at a given doping fac- 
tor, we note that an important part of the charge, about 
50% of the total hole content, is localized on the api- 
cal oxygen and lanthanum atoms. Consequently, atoms 
outside the Cu02 are also affected by doping. It was 
found in a previous work^^ on YBaCuO crystals that the 
c-axis length is related to Tc and hence to the doping 
level. The major cause of c-axis variation with respect 
to hole doping being the apical oxygen positions, our re- 
sults confirm the importance of atoms outside the plane 
in holes distribution and, indirectly, to superconductiv- 
ity. The charge on lanthanum is justified by the fact that 
even for real substitution of lanthanum by strontium, a 
little portion of the hole introduced into the structure 
should stay on the dopant. Moreover, we can see that 
the fraction of holes centered around the copper atoms is 
of the same magnitude or lower than the others atoms. 
In fact, the Bader analysis reveals that the majority of 
the holes will be injected in oxygen orbit als. This is con- 
sistent with experiments^^ and the electronic properties 
of charge transfer insulators where the oxygen orbitals lie 
just below the Fermi level. Finally, we note that hole for- 
mation increases more rapidly on in-plane oxygen than 
on other atoms. This means that for large values of dop- 
ing, hole content should increase faster in the copper 
oxide planes then in the rest of the system. The hole 
formation in Cu02 planes also explains the decrease of 
magnetic moments on the copper atoms with increased 
doping. 
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FIG. 6. Location of holes in the antiferromagnetic phase of 
doped LaCuO for [7=4 eV. Lines are guides to the eyes. 

We also conducted a Bader analysis on spin densi- 
ties for the same doping as the previous charge analy- 
sis. From this, we found that apical oxygen atoms are 
slightly spin polarized. In fact, the non-doped crystal 
shows a magnetic moment of 0.01 fiB on the apical oxy- 
gen atoms. This small spin polarization is also found in 
the 3^22 orbitals of the copper atoms. This value dou- 



bles for X = 0.25 system. The increase follows the hole 
content in these orbitals. This phenomenon, in oppo- 
sition to the antiferromagnetism in the Cu02 plane, is 
due to the direct exchange between 2pz orbitals of api- 
cal oxygen atoms and copper 'idz^ orbitals. A Wannier 
functions decomposition of the band structure for these 
systems reveals that 2pz and 'idz'2 orbitals present the 
second largest overlap. The first one being the overlap of 
the ?>dx'2j^y2 and 2px of the in-plane oxygen atoms. In ad- 
dition, we found that this polarization is noticeable only 
for values of U higher than 2 eV. 

For the stripe magnetic orders, the Bader analysis 
shows the break of translation symmetry. For the BC and 
the AC stripes, the charge modulation shows a commen- 
surability of 4 unit cells as observed in experiments.'SFor 
both types of stripes, we found that the charge modula- 
tion is practically all on the oxygen atoms. The values for 
the charges obtained for the in-plane oxygen atoms are 
reported in Fig. |2] and Fig.jS] Considering the BC stripe 
first, we found that it does not show a noticeable modu- 
lation of copper charges. The atoms in the antiferromag- 
netic parts have a charge of 27.94 electrons and copper 
atoms in the ferromagnetic domains have a smaller value 
by only 0.01 electron. Contrary to the antiferromagnetic 
phase, stripes present different atomic charges for both 
in-plane oxygen atoms. The Cu02 planes now have two 
different types of oxygen atoms. One is located between 
the copper atoms along the 8 units supercell stripe (Oi) 
and the other is parallel to the magnetic configuration 
(O2) (see Fig.jl]). This is another evidence of the symme- 
try breaking induced by stripe ordered phases. Oi atoms 
have their maximum charges between copper atoms in 
the middle of the antiferromagnetic regions. On the con- 
trary, their minimal charge, which is at the same time 
their maximum in hole content, is centered on the oxy- 
gen atoms between copper atoms in the ferromagnetic 
wall. Consequently, holes are favored in the ferromag- 
netic domains. For O2 atoms, the charge maximum and 
minimum follow the trend of the magnetization of the 
copper atoms, meaning that the larger charge values are 
found in antiferromagnetic chains with largest value of 
magnetization, which is essentially the same values as for 
the non-doped system, and the lowest charge values are 
found in the lower magnetization chains. Also, the total 
charge of O2, compared to Oi, is greater by 0.08 elec- 
tron over the 8 units cell. It indicates that BC stripes 
favor hole formation on oxygen atoms along the mag- 
netic modulation. With these results, we see that the 
antiferromagnetic region is less doped compared to the 
ferromagnetic part. A charge modulation on apical oxy- 
gen atoms is also obtained. It is due to modulation of 
magnetization on copper atoms and also to a small tilt 
of the octahedra of the Cu02 planes. 

The same analysis can be done for the AC stripe. First, 
the copper charges are almost the same for every atoms 
in the antiferromagnetic regions. For this stripe, all cop- 
per atoms have the same charge value of 27.94 electrons. 
This time, the modulation of copper magnetization is not 
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followed by a charge variation. Nonetheless, the magnetic 
ordering still splits the charge symmetry for Oi and O2. 
For Oi and O2, the modulation of charges is only sub- 
stantial for oxygens near the non magnetic walls. Partic- 
ularly, the minimum of O2 charge is on the atom bonded 
to the non magnetic copper and the minimal charges of 
Oi are found on each side of the non magnetic walls. The 
maximal charge variations are 0.01 and 0.05 electron for 
the first and second type of in-plane oxygen. In addition, 
the difference of the total charge between Oi and O2 is 
0.1 electrons over the whole stripe which is slightly larger 
than for the BC stripe. For the AC stripe, the Bader 
analysis shows that oxygen atoms around non magnetic 
copper can be compared to those of highly doped sys- 
tems, especially for O2. Finally, the difference between 
BC and AC stripes, in the Cu02 plane, is the range where 
holes are located. In fact, for AC stripes, an appreciable 
part is localized on oxygen atoms in the non magnetic 
chains whereas for BC stripes, the charge is more evenly 
distributed. 

If we now consider the bond length between copper 
and Oi atoms {x axis), we see that it varies depending 
on the charges on the neighboring oxygen atoms. For 
the BC stripe, the distance between copper and oxy- 
gen atoms along the x direction ranges between 1.89 and 
1.91 A. The maximal bond length is found for copper 
atoms linked to an Oi atom that has a charge of 9.19 elec- 
trons and to an O 2 atom with a charge of 9.21 electrons. 
The charge disparity between the two types of oxygen 
atoms surrounding the same copper atom modifies the 
bond length. The longer bond length of the Oi atom in- 
dicates a weaker bonding to the copper atom compared 
the O2 atom. Moreover, the other Oi atoms in the anti- 
ferromagnetic regions show bond lengths of 1.91 A with 
their nearest copper atoms, the difference in charge be- 
tween Oi and O2 atoms being insufficient to noticeably 
modify the bond lengths. On the contrary, the 1.89 A 
bond length found in the ferromagnetic domains are as- 
sociated to Oi atoms having the lowest atomic charge. 
The reason is that, due to the symmetry in the y axis, 
CU-O2 bond length are fixed to 1.9 A as in the antifer- 
romagnetic zone. The O2 atom having a smaller atomic 
charge in ferromagnetic regions than in antiferromagnetic 
regions, the fixed bond length of O2 atoms prevents this 
atom from being closer to the copper atom. As a re- 
sult, it favors the reduction of the bond length between 
copper atoms and its closest Oi atom. Finally, the BC 
stripe phase affects simultaneously the copper magnetic 
moments, the charge distribution and the Cu02 plane 
structure to lower the total energy of the system. In the 
same way, bond lengths in the AC stripe are modulated 
along the x direction. For this phase, the maximum bond 
length is still 1.91 A and is found in the non-magnetic 
wall. As discussed before, the larger bond length ob- 
tained for Cu-Oi link is explained by difference between 
Oi and O2 atomic charges, when the O2 atom has the 
maximal atomic charge (9.21 electrons). In addition, the 
smallest bond length is also found in the non-magnetic 



wall. In fact, the AC stripe presents the lowest atomic 
charge for in plane oxygen atoms (9.16 electrons) when 
compared to the BC stripe. This important concentra- 
tion of hole localized on an O2 atom, combined to the 
fixed CU-O2 bond length (1.9 A) result in a bond length 
of 1.88 A, similarly to the BC stripe. The smallest bond 
length found in the AC stripe, compared to its analogue 
in the BC stripe, is accounted by the higher hole con- 
tent in the O2 atoms of the AC stripe. The fact that 
BC stripe is more stable than the AC stripe indicates 
that the structure favors shorter bond lengths modula- 
tion along the x axis. The proximity of the Oi atom and 
the copper atom in the 1.88 A bond of BC stripe come 
with a small increase in the total energy of the system 
due to Coulomb repulsion, which explains, partially, the 
higher energy of this phase compared to the BC stripe. 

The Bader analysis on the BC structure can be com- 
pared to the antiferromagnetic state of the non-doped 
system to give a more precise location of holes in this 
stripe. The study of Bader charges in antiferromagnetic 
phases showed that copper atoms only hold a fraction 
of the total holes injected in the structure. This is still 
the case for BC stripes. In fact, copper atoms possess a 
slightly larger hole fraction, 0.18 hole, compared to the 
antiferromagnetic phase with the same doping. It can 
be accounted by the variation of the copper magnetic 
moments, which tend to result in a increased hole con- 
centration on low magnetization atoms. In the case of 
in-plane oxygen atoms, they hold a major part of the to- 
tal hole content, 0.37 hole, and it is distributed almost 
evenly with a larger hole content on Oi oxygen atoms 
between the two copper atoms in the ferromagnetic re- 
gion. In this type of stripe, holes are more localized on Oi 
oxygen atoms as compared to O2 oxygen atoms by about 
33%. As mentioned before, the modulation of magneti- 
zation induces a reorganization of holes in the system. 
The charge analysis revealed that the BC stripe creates 
a charge modulation for both types of in-plane oxygen. 
Moreover, this charge modulation creates hole rich and 
hole poor regions in the CUO2 plane. 

The AC stripe has a smaller amount of hole on their 
copper atoms, compared to BC configuration, but a 
larger portion of hole is located on oxygen atoms in the 
Cu02 planes. The portion of hole doping on copper 
atoms is reduced to 0.14, comparable to the antiferro- 
magnetic phase. On the contrary, the in-plane oxygen 
atoms share 0.43 hole, which is higher than in the BC and 
antiferromagnetic phases. The high value of hole content 
in the non magnetic wall and oxygen atoms nearby cre- 
ates a large concentration of hole in this region. In fact, 
an amount of 0.15 hole is found on these oxygen atoms. 
This represents a quarter of the total hole contained in 
the Cu02 plane. The localization of hole is more impor- 
tant in this stripe, where non magnetic walls hold a large 
fraction of hole. 

Finally, for x=0.125, the BC stripe order is found to 
be the ground state of the system. From this and the 
charge analysis, we can conclude that this doped phase of 
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La2Cu04 is more stable than for configurations inducing 
hole localization in small regions of the Cu02 plane. 



IV. CONCLUSION 

To conclude, several configurations of La2Cu04 were 
studied, for hole doping ranging from to 0.25, using 
DFT+U. In the cases where supercells allow a charge 
order half the length of the magnetic order, a bond cen- 
tered stripe is found to be the ground state of the sys- 
tem. The cases x = \ and ^ = | satisfy this condition. 
For a doping of ^, no stripes are energetically favored 
and the lowest energy system is the antiferromagnetic 
state. This configuration does not allow the same pat- 
tern as X = \ OY X = |. Moreover, a charge investi- 
gation on doped phases of antiferromagnetic La2Cu04 
indicates that holes are located mainly on oxygen atoms. 
Also, study of the spin density found that apical oxy- 
gen atoms have a small magnetization, caused by the 
overlap of 2pz and 'idz'^. Furthermore, hole formation is 



found to be greater for in-plane oxygens, compared to 
copper atoms. For stripes, the modulation of magnetic 
moments on the copper atoms can be understood by the 
charge reorganization of the oxygen atoms of the Cu02 
planes. This means that stripe ordered phases deform 
lightly the crystal to favor a different magnetic order than 
the usual antiferromagnetic conformation. Also, it was 
shown that stripe orders favor the vicinity of hole poor 
and hole rich regions. This reorganization of the elec- 
tronic density should induce modification to the phonon 
spectra for doping factor allowing stripe phases. 
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